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Abstract 

The model-independent piecewise parametrizations (0-spline, linear-spline and cubic-spline) are 
used to estimate constraints of equation of state of dark energy (wd e ) from current observational 
data ( including SNIa, BAO and Hubble parameter ) and the simulated future data. A combination 
of fitting results of Wd e from these three spline methods reveal essential properties of real equation 
of state Wde- It is shown that Wde beyond redshift z ~ 0.5 is poorly constrained from current 
data, and the mock future ~ 2300 supernovae data give poor constraints of Wd e beyond z ~ 1. 
The fitting results also indicate that there might exist a rapid transition of Wd e around z ~ 0.5. 
The difference between three spline methods in reconstructing and constraining Wd e has also been 
discussed. 
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I. INTRODUCTION 



The current expansion of the universe is found to be accelerating jl, 2|, one of the possible 
explanations for this is the existence of dark energy (DE), whose energy density is dominant 
in the universe and its present equation of state Wd e o is less than —1/3. At present, it is 
still fair to say that one knows a little about the nature of DE. From current astronomical 
observational data, we can obtain some properties of the equation of state Wd e (the ratio 
of the pressure and energy density of DE). To fit the observational data, one has to first 
assume a form of Wd e - Some forms of Wd e have been used in the literature. For example, the 
CPL parametrization Wfe(z) = w + w a z/(l + z)^, 4] and the redshift expansion, Wfe{z) = 
Wo + w zZ [s 8 ] . It has been found that Wd e is very close to —1 and should be varying very 
slowly with redshift (if any). In the most cases, the cosmological constant with w = — 1 is 
still favored within 2a confidence level (C.L.). Of course the fitting results are dependent on 
the parametrization forms adopted. Usually a parametrization form is only suitable to mimic 
one type of Wd e {z). For example, the CPL parametrization can describe linear and smooth 
Wde well but is hard to reconstruct Wd e with oscillations or rapid transitions. On the other 
hand, several model-independent methods have also been proposed to reconstruct Wd e (z) |9j- 



131 ]. e.g., the uncorrelated band-power estimate 



linear-spline method 



16r|l8|, wavelet approach 



i.e., 0-spline) [14j, cubic-spline method 



191 ]. and Gaussian process modeling 




etc. Most of mo del- independent methods have a piecewise Wd e (z): 

Wde(z) = Wi(z), Zi< Z < Zi+i, (1) 

where Wi(z) is a simple function of redshift z, e.g., Wi(z) is just a constant in each bin for 
the 0-spline method. 

Different mo del- independent methods should give different but consistent fitting results of 
DE. In this paper, we would like to get constraints of Wd e from present and next generational 
observations by using three piecewise parametrizations: 0-spline, linear-spline and cubic- 
spline. The difference among the three spline methods in constraining Wd e will be analyzed. 
It is shown that each spline method is only suited to certain types of Wd e - A combination of 
constraints of Wd e from the three spline methods should help to get real properties of Wd e - 

At first, the constraints of Wd e from current observational data will be obtained. It is 
found that the present constraints on Wde are very weak beyond z ~ 0.5, because in the 
higher redshift region there are less data points and the effect of DE on the expansion of the 



universe is weaker. Moreover, we see from the fitting results from 0-spline and linear-spline 
methods that there might exist a rapid transition of Wde around z — 0.5. To estimate the 
constraints on Wd e from next generation observations, ~ 2300 SN data as forecasted for a 
space mission like SNAP/JDEM 23|, |24| are simulated. The mock data are simulated from 
two fiducial models: one with a smooth Wd e and the other with a rapid transition u^ e around 
z = 0.5. In this case Wd e is poorly constrained after z ~ 1. 



II. THE CONSTRAINTS ON w de FROM PRESENT DATA 



We use three spline methods to get constraints of Wd e from present observational data, 
which include type la supernovae data (SNIa), baryon acoustic oscillation (BAO) and ob- 
servational Hubble data. The best-fitting parameters and their errors will be obtained by 
using the Markov Chain Monte Carlo (MCMC) method. In general, the errors of Wd f . of 
different bins are correlated. We will adopt the decorrelated method proposed in 10j to 
obtain errors of a new parameter Q(z). The new parameter Q(z) is defined by transforming 
the covariance matrix of Wd e , so that the errors of Q(z) are uncorrelated and do not entangle 
in each bin. 



A. Observational data 

Since it is hard to get strong constraints of Wd e in high redshift region, we will focus 
on constraining Wd e in the region z G [0, 0.9] only in this work. The data sets we adopt 



are SNIa Union2 data 
from 



28 



25], BAO distance parameter A 



261 ] and observational Hubble data 



30j. Those data points with z > 0.9 in these data sets will be abandoned. As a 



result, in our calculations, only 519 SNIa data points, 11 Hubble data points and one BAO 
data point are used. In this case, one needs not to assume the form of Wd e {z > 0.9) and to 
consider correlations between Wd e (z > 0.9) and Wde(0 < z < 0.9). 
The cosmological parameters are fitted with the SNIa data 31] by 



XsN 



519 

£ 

i=l 



07 



(2) 



where the theoretical distance modulus 



ji th (z) = 5 log 10 Dl + fMo 



(3) 



For a flat Friedmann-Robertson- Walker universe, the luminosity distance is 

D L = (1 + z) [ dx/E(x) , (4) 



where 

E 2 (z) = n m0 (i + zf + (1 - a m0 )F(z). (5) 

Here f2 m0 is the current fractional matter density of the universe and function F(z) depends 
on the parametrization of w^z): 

F(z) = e 3 X 1 -^ d *. (6) 
One can expand Eq.(j2J) with respect to /x as 

XSN = a + 2 Vo + Wo ( 7 ) 

where 

_ \lHh{zj;no = 0) ~/4] 2 
i=i ' 

519 

6 = ^ /i^(^;/^o = 0)-^ fe ^ 
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C 



The Xsn nas a minimum with respect to 

xl N = a- b 2 /c . (9) 
This way the nuisance parameter /i is reduced, in this work we will adopt Xsn instead of 



Xsn- 



The BAO distance parameter A is the measurement of BAO peak in the distribution of 



SDSS luminous red galaxies [26 1 



A = ^E-^)[^f-^f. (10) 



The value of A is determined to be 0.469(n s /0.98 ) 0,35 ±0.017, where n s = 0.963 is the scalar 
spectral index, which has been updated from the WMAP7 data [27]. The x\ i s defined as 



(A - 0.472)' 



Xa _ 0172 • (11) 
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z 





0.1 


0.17 


0.27 


0.4 


0.48 


0.88 


0.9 


0.24 


0.34 


0.43 


h 


0.738 


0.69 


0.83 


0.77 


0.95 


0.97 


0.9 


1.17 


0.7969 


0.838 


0.8645 




0.024 


0.12 


0.08 


0.14 


0.17 


0.6 


0.4 


0.23 


0.0232 


0.0296 


0.0327 



TABLE I. 11 observational Hubble data with their redshifts from 



28|43C] 



The observational Hubble data can be obtained by using the differential ages of passively 
evolving galaxies 

1 Az 



H 



1 + z At 



(12) 



Here we adopt 11 observational Hubble data points with z < 0.9 from 
points are summarized in Table [B The Xhub * s defined as: 



28 



-|30|. Those data 



2 _ \hth{zi) - h ob (zi)] 2 
Xhub — 

i=i 



(13) 



where h = if/100 km ■ s 1 • Mpc 1 . 

Finally the total Xt f° r three kinds of observational data is the sum of them: 



2 _ ~2 ,2,2 

Xt — Xsn ' Xa ' Xhub- 



(14) 



B. Methodology and Results 

Now by the piecewise parametrization approaches for Wd e , we get the constraints of Wd e 
from the observational data mentioned in the above. We will use the 0-spline, liner spline 
and cubic spline method, respectively. The best-fitting w de (z), Q(z) and their la, 2a C.L. 
errors obtained by using three spline methods are plotted in Fig. [1] . One can see that the 
constraints of Wd e from the three methods are consistent with each other. It is shown that 
the constraints of Wd e beyond z ~ 0.5 are much weaker than those within z ~ 0.5, and 
Wde = — 1 (the cosmological constant) is still consistent with present constraints of Wde at 
2a C.L. In addition, there might exist a rapid transition of Wd e around z = 0.5, which is 
particularly evident in the fitting results from the 0-spline and linear-spline methods. 
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FIG. 1. The best fitting Wd e (z), Q(z) and their la and 2a C.L. errors from the current observational 
data (SN+A+HUB). The black dotted curves stand for w = — 1. The three figures in left panel are 
for the correlated Wd e (z) and the three figures in right panel are for the uncorrelated parameters 
Q(z). The upper two figures are from the 0-spline method, the middle two from the linear-spline 
method and the bottom two from the cubic-spline method. Note that some best fitting results on 
Wde and their errors are not included in figures when they are beyond the range of figures. 
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1. 0- spline method 



In this method, one divides the redshift region under consideration into n bins, and sets 



where Wi are constants and Z{ are divided positions of bins. Since the number and precision 
of present data are still not sufficient enough, at first we divide the redshift region [0, 0.9] 
into only two bins and treat the divided position z\ as a free parameter of the model. By 
fitting the model with data, we find the best value of divided position z\ = 0.62. Thus we 
have 



In the fitting process we have assumed a prior —20 < w%. Otherwise W2 would go to a 
very large minus value in MCMC procedure and thus the downward error of w 2 would be 
extremely large. Of course the prior will not qualitatively affect the results and conclusions 
of this paper. It is found that there is almost no difference between the errors of u>; and 
corresponding Q t . This shows that the correlation between W\ and W2 is extremely small in 
this case. 

The fitting results are shown in Table [TT] and Fig. [TJ It indicates that the constraints of 
Wde are very good in the whole first bin, since here we have assumed Wd e to be a constant 
and the errors of Wd e are averaged in each bin. For the second bin, the errors of W2 are very 
large, particularly its downward error, as expected. Of course the fitting results depend on 
the number of bins and the divided manners. Here the width of the second bin is relatively 
small and there are only 103 data points in this bin, but the main reason for the weak 
constraint of w 2 (wd e in the second bin) is due to the high redshift; in the lower redshift 
region, the same width of bin and the same number of data points can give much better 
constraints of Wde- 

Next we divide the readshift region z G [0, 0.9] into three bins to see whether there exist 
more structures of Wd e (z) in this region. Two divided positions are also treated as free 
parameters, the positions of them from the best fitting are: z\ = 0.45 and z 2 = 0.50. In this 



W de (Zi <z< z i+1 ) = w. 



(15) 




(16) 
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h 


^m0 


w\ and Q\ 


u> 2 and Q 2 


y 2 ■ 


fl 799+0.012+0.025 
u - ' zz -0.014-0.028 


n 97C+0. 022+0.042 
u - z ' d -0. 015-0.033 


-, nA +0.06+0. 12 
1 - u ^-0.06-0.14 

1 r.A+0. 06+0.12 


. c-9+3. 76+4.40 
^• oz -15. 38-15. 47 

a a 9+3.42+4.04 
^•^-15.00-15.05 


499.832 



TABLE II. The best-fitting values with la and 2a C.L. errors of parameters from present data in 
the case with the 0-spline method. 



case, the second bin is relatively narrow. The best-fitting Wd e (z) is found to be: 



W d e{z) 



-1.06 , < z < 0.45 
5.82 , 0.45 < z < 0.5 
-5.40 , 0.5 < z < 0.9 



(17) 



This results are plotted in Fig. CEJ but the result for the second bin u? 2 =5.82 is not included 
because it deviates far from Wd e in other bins. The appearance of the narrow second bin 
and the large deviation of w-i from the values in other bins might imply that there is a rapid 
transition of Wde around z ~ 0.5. Mock future data with a rapid transition Wde around 
z ~ 0.5 is simulated in the next section and a similar fitting result from the data is found. 



2. Linear-spline method 

In this case we set Wde(z) as 

w(zi) — w(zi-x) , 

w de {Zi-i < z < Zi) = w{zi_i) H {z-Zi-t), (1 < i < n) (18) 

Z% Z{ „ ] 

where n is the number of bins and zq = 0. Still, the redshift region z G [0, 0.9] is divided 
into two bins and the best fitting result gives the divided position to be Z\ = 0.56, which 
is very close to the one given by 0-spline method. Here the prior —20 < u>(0.9) has been 
assumed. The best fitting values of the parameters w(0), w(0.56) and w(0.9) and their 
errors are shown in Fig. [1] and Table III. It is shown that the errors of Wd e in the second bin 
(especially the downward error) increase quickly with redshift z. At z — 0.9, the constraint 
of Wde becomes extremely weak. 

The best-fitting Wd e {z) in [0, 0.9] with three bins is also shown in Fig. [U The best fitting 
results for w de are: w(0) = -0.78, w(0.42) = -1.68, iu(0.43) = 4.28 and w{0.9) = -20. 
Note that here we have also assumed the prior —20 < tu(0.9). One can see from the figure 
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that there is also a rapid transition of Wd e around z ~ 0.5 in this best fitting, as the case of 
the 0-spline method. 



h 


^m0 


w(0) and Q(0) 


w(0.56) and Q(0.56) 


w(0.9) and Q(0.9) 


x 2 ■ 

^mvn 


n 79-1+0.012+0.026 
u - '^ i -0.013-0.027 


n 97fi+0.026+0.046 
U.Z(U_ Q .015-0.035 


n qr+0.16+0.31 
u - yo -0.15-0.31 

i m +0.06+0.13 
1 - ui -0.07-0.15 


i 9A +0.37+0.75 

-L-^-o^-o.sg 

i -, o+0. 12+0.25 
1 - lo -0. 19-0.38 


Qo+6. 13+7.85 
^• yo -15.06-15.06 

A /|f:+4.97+6.69 
^•^0-13.02-13.08 


499.692 



TABLE III. The best-fitting values with \a and 2a C.L. errors of parameters from present data in 
the case with the linear-spline method. 



3. Cubic-spline method 

To use the cubic-spline method, we divide z G [0, 0.9] into three bins with the fixed 
divided positions as: Z\ = 0.3, z 2 = 0.6. The results are shown in Fig. [1] and Table IV. In 
this case, once again, the constraints of Wd e in the last bin are very weak and the errors of 
Wde are extremely large. The errors of the uncorrelated parameter Q(0.9) are much smaller 
than those in other two methods but the errors of other Q's are larger than those in other 
two methods, because in the cubic-spline method Wd e in different bins are highly correlated. 



h 


f^mO 


w(0) and Q(0) 


to(0.3) and Q(0.3) 


w(0.6) and Q(0.6) 


u;(0.9) and Q(0.9) 


Xmin 


n 7on+°- 013 + - 026 

U. 1 4U_ ,013-0 .027 


n 970+0. 026+0. 049 
u - z '°-0. 015-0. 035 


n /in+°- 42 + - 91 

u - w -0.63-1.17 

n qq+O.12+0.25 
u - JZ -0. 15-0.28 


1 nc+0.16+0.32 
l.Ud_ .20-0.43 

i n8+°- 06 + - 14 

J - UO -0.08-0.17 


1 cfi+0.92+1.60 
-L- 00 -!. 10-2. 55 

_n QC+0.58+X.24 
U.»O_ .40-O.8O 


Q oo+7. 85+12. 66 
y - oo -10. 09-10. 11 

9 00 + 1.17+2.59 
z -°°-1.94-2.10 


498.942 



TABLE IV. The best-fitting values with la and 2a errors of parameters from present data in the 
case with the cubic-spline method. 



III. CONSTRAINTS OF w de FROM FUTURE DATA 

To see the constraint ability on Wd e from future observational data, we adopt the charac- 



teristics of a SNAP-like JDEM survey 23) to simulate the future SN data with 0.1 < z < 1.7, 



which include 1998 SN data points. The redshift distribution of the mock SN data is shown 



z— s> 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.1 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


Nbin 


300 


35 


64 


95 


124 


150 


171 


183 


179 


170 


155 


142 


130 


119 


107 


94 


80 



TABLE V. The redshift distribution of 1998 SN data with 0.1 < z < 1.7 for a SNAP-like JDEM 



survey 



23J and 300 SN data with z < 0.1 from the NSNF 



24|. The redshifts in the first row are 



the upper limits of each bin. 



19 



in Table V, in which 300 supernovae with z < 0.1 
bin as shown in Table V, SN's are assumed to be uniformly distributed 

We will use two fiducial models to simulate the mock data: 
Model I: one assumes a slowly varying equation of state for DE: 

300000 



241 ] are also included. In each redshift 



Wde(z) 



-0.8 - 



e 22/(i+*) + 600000 ' 

Model II: one has the equation of state with a rapid transition around z ~ 0.5 : 

3 x 10 14 



-0.8 



(19) 



(20) 



e 100/(l+z)-30 _|_ 6 x 10 14 ' 

In Fig. |2] two fiducial Wd e (z) are plotted. The form of Wd e {z) in these two models is also 
adopted in (20]. Both fiducial models have h = 0.72 and Q m o = 0.28. Now the distance 



modulus of ~ 2300 SN can be simulated and the corresponding errors are assumed as 



19 



23| : 



a 



aiz 



obs 



+ dm 2 , 



bin 



(21) 



where a a b s = 0.15, dm = 0.02z/z max and z max is the maximum redshift (here z max = 1.7). 
To simulate the effect of other future observations and alleviate the degeneracy between 
fl m o and Wde, we add a prior Q m o = 0.28 ± 0.03. In all three spline methods, the redshift 
region (0, 1.7) will be divided into three bins. Two divided positions are still treated as free 
parameters in the cases of 0-spline and linear-spline methods, while the divided positions of 
bins are fixed by hand in the case of cubic-spline method. All fitting results are shown in 
Fig. |2j One can see that three spline methods give consistent results and the future mock 
data give poor constraints of Wd e beyond z ~ 1. 



A. 0-spline method 



For the fiducial model I, the region (0,1.7) is divided into three relatively uniform bins in 
the best fitting model, with the divided positions: Z\ = 0.55 and z 2 = 1.27. The errors of 

10 




FIG. 2. The best fitting Wd e (z) and their la and 2a C.L. errors from the mock future data. The 
dashed curves are for the fiducial models, and the black dotted curves stand for the cosmological 
constant w = —1. The three figures in left panel are for model I and the three figures in right 
panel for model II. The upper two figures are from the 0-spline method, the middle ones from the 
linear-spline method and the bottom two figures from the cubic-spline method. 



Wde increase rapidly with redshift. In the last bin, the best-fitting Wd e deviates far from the 
fiducial model and the errors of Wd e are relatively large. 

For the rapid transition model II, the best fitting divided positions are Z\ = 0.52 and 
Z2 = 0.53. In this case the second bin is extremely narrow because of the rapid transition of 

11 



the fiducial Wde around z ~ 0.5. The best fitting Wde in the second bin is w 2 = 0.67 which 
deviates from the fiducial Wde beyond 2a C.L. . The errors of Wd e in the last bin are much 
better than those for model I, since here the width of last bin is much larger than that in 
model I. 

B. Linear-spline method 

In this case the situation is very similar to the case of the 0-spline method. For model I, 
the best fitting divided positions of bins are Z\ = 0.3 and z 2 = 1.35. It is shown that the 
fiducial Wde can be well reconstructed until z ^ 1. For model II, the best fitting divided 
positions are Z\ = 0.43 and z 2 = 0.62. Here the constraints of Wd e at high redshift are much 
better than those for model I, since the last bin here is much larger than that for model I. 
It can be seen that the rapid transition of Wd e in model II can be well reconstructed by the 
linear-spline method, as shown in Fig. [21 though the width of second bin is very narrow. 
The reconstructed Wd e and its errors here have finer structures than those in the 0-spline 
method. 

C. Cubic-spline method 

In this case we fix the divided positions of three bins as z\ = 0.5 and z 2 = 1.0. The 
errors of Wd e still increase rapidly with redshift. The fitting results for model I and model II 
indicate that this method can reconstruct the slowly varying Wd e well, but it is not good in 
reconstructing the equation of state with rapid transition. For model I the errors of Wd e are 
consistent with those from other two methods, but for model II the errors are much larger 
than those from other two spline methods. 

D. Result analysis 

For model I, it is shown that Wd e can be well reconstructed up to z ~ 1. For model II, 
since the width of the last bin are always much larger than that in model I, the errors of 
Wde in the region beyond z ~ 1 are much smaller than those in model I (except for the case 
of the cubic-spline method, which is not good at describing a rapid transition Wd e )- This 
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means that the fitting results depend on the divided manner of redshift bins. In our case, 
we treat the divided positions as free parameters and then fix their values to the best-fitting 
values (in 0-spline and linear-spline methods) or just divide the redshift region uniformly (in 
the cubic-spline method). In this case, the width of the last bin is always not large enough 
to get strong constraints of Wd e in that bin. One may use other ways to divide redshift, even 
setting the last bin large enough by hand. But a large bin usually will lead to a lose of fine 
structure of Wd e in this case. 

In the case of the 0-spline method, Wd e and its errors are averaged in each bin, while 
in the cases with other two spline methods, the errors of Wd e increase with redshift inside 
each bin and the reconstructed Wd e always have finer structures than that from the 0-spline 
method. For the cubic-spline method it is good at reconstructing the slowly varying Wde, 
but not the case with rapid transitions, and the errors of Wd e are highly correlated. 

IV. CONCLUSIONS 

We have studied the constraint ability on the equation of state Wd e of dark energy from 
the present and simulated future observational data by piecewise parametrization with the 
0-spline, linear-spline and cubic-spline methods, respectively. Three spline methods give 
consistent results of Wd e - 1) the cosmo logical constant Wd e = — 1 is still consistent with 
present data at 2a C.L.; 2) current data can constrain Wd e well up to z ~ 0.5 and the future 
(mock ~2300 SN) data can constrain u> de well up to z ~ 1; 3) in high redshift region, the 
downward errors of Wde are always much larger than the upper ones; 4) the fitting results 
from current data by using the 0-spline and linear-spline methods indicate that there might 
exist a rapid transition of Wd e around z ~ 0.5. 

There are also differences among the fitting results from the three spline methods. With 
the 0-spline method Wd e and its errors get averaged in each bin, and thus it always gives 
poor structure of Wd e - Therefore this method is suited to be used to confirm whether Wd e 
is a constant (including -1) or not. The linear-spline and cubic-spline methods give finer 
structure of Wd e than the 0-spline method. The linear-spline method can reconstruct almost 
all types of Wd e in principle, but the reconstructed Wd e is always not smooth at the divided 
positions of bins, which will lead to deviations of Wd e from the real Wd e around the divided 
positions. Thus the linear-spline method is suited to reconstruct non-smooth Wd e (z), and the 
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positions where Wd e suddenly changes can be accurately determined. For the cubic-spline 
method, one needs not to search for the best-fitting divided positions of bins but the redshift 
region must be divided into at least 3 bins. It is shown that the cubic-spline method is not 
good at reconstructing Wd e with rapid transitions and the errors of Wd e at different bins are 
highly correlated. The cubic-spline method is therefore suited to reconstruct a smoothly 
varying Wd e - Basically, a combination of the fitting results from the three spline methods 
can reveal the real Wd e - 

The fitting results are also affected by divided manners of redshift bin. Usually a larger 
width of one bin will lead to a stronger constraint of Wde there, but fine structure of Wd e will 
be lost. At present, the number and precision of observational data are still not sufficient 
to obtain both strong constraints and fine structures of Wd. P .. In particular, to constrain Wde 
at high redshift, a large number of data will be required Q. 
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